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Three-Dimensional  Simulations  of  a  Gas/Gas,  Hydrogen/Oxygen  Engine 

Mark  Archambault 

Air  Force  Research  Laboratory 
Propulsion  Directorate 
Space  and  Missile  Propulsion  Division 
Edwards  AFB,  CA  93524 

Oshin  Peroomian 
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Westlake  Village,  C A  91361 

The  overall  objective  of  this  research  is  to  establish  a  design  methodology  for  gas/gas  injectors.  This  paper, 
however,  focuses  on  a  computational  methodology  to  efficiently,  accurately,  and  robustly  obtain  high- 
fidelity  solutions  of  combusting  rocket  engine  flows  to  gain  a  knowledge  and  understanding  of  their 
features.  To  that  end,  simulations  of  a  single-element,  shear-coaxial,  H2/O2  engine  are  being  performed  to 
characterize  its  flowfield  and  to  validate  the  CFD++  flow  solver  for  this  class  of  problems.  Previous  work 
has  focused  on  obtaining  solutions  on  a  grid  three  to  four  times  finer  than  those  reported  by  other 
researchers  and  resolving  numerical  issues  that  reduce  the  computational  efficiency  of  this  inherently 
unsteady  flow.1-4  Comparisons  of  steady  and  averaged  time-accurate  solutions  have  also  shown  that  a 
steady  solution  may  not  provide  an  accurate  depiction  of  the  combusting  flow  field  over  time  (Fig.  1). 

Other  simulations  have  shown5  that  flow  features  unique  to  an  experimental  configuration,  such  as  a 
nitrogen  curtain  purge  used  to  cool  the  optical  access,  can  influence  both  the  experimental  and 
computational  results.  Figure  2  shows  that  when  the  nitrogen  curtain  purge  present  in  the  experiment  is 
modeled,  the  predicted  hydrogen  profile  more  closely  matches  the  data.  This  is  due  to  the  fact  that  when 
the  nitrogen  is  present,  the  hydrogen  is  unable  to  radially  diffuse  to  the  engine  walls  as  quickly  as  when  the 
nitrogen  is  absent.  It  is  clear  that  the  nitrogen  has  had  an  influence  on  the  experimental  data  and  both  the 
experimenter  and  the  modeler  should  take  care  when  interpreting  their  results. 

Finally,  a  comparison  of  turbulence  models  has  provided  evidence  that  the  simulation  is  sensitive  to  which 
turbulence  model  is  employed.5  Figure  3  shows  profiles  of  hydrogen  mole  fraction  and  mean  axial  velocity 
generated  by  various  turbulence  models.  It  is  clear  from  the  figures  that  the  results  are  sensitive  to  the 
selection  of  the  model.  In  this  case,  we’ve  found  that  a  linear  k-epsilon  model  best  predicts  the  data.  This 
is  likely  caused  by  higher  values  of  turbulent  viscosity  produced  by  the  linear  k-epsilon  model,  allowing  for 
increased  mixing  with  the  oxygen. 

The  current  work  continues  these  studies  by  examining  the  results  of  three-dimensional  calculations  on  this 
same  configuration.  The  first  simulation  will  determine  whether  there  exists  a  three-dimensional  flow 
structure,  such  as  a  helical  shear  layer,  when  using  the  axisymmetric  geometry.  Steady  and  time-averaged 
results  (including  the  nitrogen  curtain  purge)  will  be  discussed  and  compared  with  experimental  data. 
Additional  simulations  will  focus  on  biasing  the  GOX  post  such  that  it  is  either  non-concentric  with  the 
hydrogen  annulus,  or  non-parallel  with  the  hydrogen  stream.  The  purpose  of  these  calculations  will  be  to 
determine  if  manufacturing  errors  of  this  type  could  have  a  significant  impact  on  the  flow,  and  thus  the 
efficiency  of  the  engine.  In  addition,  this  will  also  provide  us  an  initial  step  in  investigating  whether 
intentional  GOX  post  biasing  could  be  used  as  a  means  to  control  propellant  mixing. 
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Figure  1.  OH  concentration  contours  from  (a)  instantaneous, 
(b)  averaged  time-accurate,  and  (c)  “quasi-steady”  solutions. 
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Figure  3.  Profiles  of  (a)  hydrogen  mole  fraction  and 
(b)  mean  velocity  with  various  turbulence  models. 
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